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Abstract. Power density spectra (PDS) that are characteristic of low mass X-ray binaries (LMXBs) 
have been previously reported for M3 1 X-ray sources observed by XMM-Newton. However, we 
have recently discovered that these PDS are false positives resulting from the improper manipu- 
lation of non-simultaneous lightcurves. The lightcurves produced by the XMM-Newton Science 
Analysis Software (S AS) are non-synchronised by default. This affects not only the combination of 
lightcurves from the three EPIC detectors (MOS1, MOS2 and pn), but also background subtraction 
in the same CCD. It is therefore imperative that all SAS-generated lightcurves are synchronised by 
time filtering, even if the whole observation is to be used. We combined simulated lightcurves at 
various intensities with various offsets and found that the artefact is more dependent on the offset 
than the intensity. While previous timing results from M3 1 have been proven wrong, and also the 
broken power law PDS in NGC 4559 ULX-7, XMM-Newton was able to detect aperiodic variability 
in just 3 ks of observations of NGC 5408 ULX1. Hence XMM-Newton remains a viable tool for 
analysing variability in extra-galactic X-ray sources. 



INTRODUCTION 

The variability and spectral properties of Galactic low mass X-ray binaries (LMXBs) are 
well known to depend more on the accretion rate than the nature of the primary [neutron 
star or black hole, 1]. At low accretion rates, the power density spectra (PDS) may be 
characterised by broken power laws, with the spectral index, y, changing from ~0 to ~1 
at some break frequency in the range 0.01-1 Hz [1]; such variability has a r.m.s. power 
of ~ 10-40% [e.g. Si. We describe such PDS as Type A [3]. At higher accretion rates, the 
PDS is described by a simple power law with y ~1 and r.m.s. variability <10% JI10]- 
We refer to these PDS as Type B |0]. 

Type A variability is characteristic of disc-accreting XBs, and any X-ray source that 
exhibits such variability may be identified as an XB, rather than a fore-ground star or 
background active galaxy [see e.g. 4]. Such variability has been reported by Barnard et 
al. in XMM-Newton observations of M31 [see e.g. IHSELB]. Type A PDS have also 
been reported for ultraluminous X-ray sources (ULXs) in NGC 4559 ED and NGC 5408 
§]. 
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FIGURE 1. XMM-Newton combined MOSl+MOS2+pn lightcurves of r3-61 in the 0.3-10 keV band; 
the combination of lightcurves produced by Method I and Method II are shown in the top and bottom 
panels respectively. The best fit line of constant intensity to the Method I lightcurve gives ^ 2 /dof = 
408/238, while the best fit to the Method II lightcurve gives £ 2 /dof = 280/239 



XMM-NEWTON LIGHTCURVES 



However, we have now discovered an additional, artificial origin for such stochastic 
variability, caused by the improper addition of lightcurves from the three EPIC instru- 
ments on board XMM-Newton (MOS1, MOS2 andpn), where the lightcurves are not 
synchronised. This artefact is discussed in detail in 

XMM-Newton has observed M31 extensively, and Barnard et al. have 
conducted a survey of bright sources in the central region, over four observations. For 
each observation of every source, 0.3-10 keV lightcurves were obtained from source 
and background regions in the MOS1, MOS2 and pn images, using 2.6 second binning. 

We have recently discovered that the results of manipulating lightcurves produced 
with SAS tool evselect using lcmath produces different results, depending on how the 
time selection is performed. The three methods of time selection may be summarised as 
follows 

• Method I: using the expression "(TIME in [tstart.tstop])" in evselect 

• Method II: adding the keywords TLMIN1 = t start and TLMAX1 = tstop to the 
events files for each detector, using the fparkey FTOOL 

• Method III: filtering each lightcurve with the evselect keywords "timemin =tstart 
timemax=tstop , 

Methods II and III are equivalent. Vitally, they are not equivalent to Method I, although 
this is not mentioned in any documentation. 

In Fig. d we compare MOSl+MOS2+pn lightcurves of the M31 X-ray source 
XMMUJ004207.6+411815 (r3-61) from the 2002, January 6 XMM-Newton observa- 
tion, extracted using Methods I (top) and II (bottom). The best fit line of constant inten- 
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FIGURE 2. Error vs. intensity for data points in MOSl+MOS2+pn lightcurves of r3-61 combined 
using Method I (left) and Method II (right). In each case, a curve represents Gaussian errors. Clearly, 
the lightcurves produced using Method I were incorrectly added, as the errors are mostly smaller than 
Gaussian; meanwhile the Method II lightcurve has Gaussian errors, showing that it has been correctly 
treated. 



sity to the Method I lightcurve yields ^ 2 /dof=408/239, while the best fit to the Method 
II lightcurve gives ^ 2 /dof=280/240. The best fit intensity is 0.102 count s _1 for both 
lightcurves. Consequently, we analysed each lightcurve to determine which, if any, was 
correct. 

We first examined the MOS1, MOS2 and pn source lightcurves individually before 
background subtraction, testing whether the errors were simply Gaussian. We binned 
each lightcurve to 100 s and obtained the count rate and error for each bin of the result- 
ing lightcurve. We found that each of the input lightcurves had been assigned Gaussian 
uncertainties (uncertainty o = N^ 5 for N counts in a time bin), so the correct combined 
lightcurve would have Gaussian uncertainties also. In fact, Gaussian uncertainties un- 
derestimate the confidence limits at low count rates lfl2ll . 

The rate vs. errors for the combined MOSl+MOS2+pn lightcurves using Method I 
and Method II are presented in the left and right panels respectively of Fig [2j It is clear 
that the Method II lightcurve has been correctly produced, as the errors are Gaussian; 
however, the errors for the Method I lightcurve are unfeasibly small, as they are smaller 
than the Gaussian limit; as a result, the Method I lightcurve is artificially variable. 

Closer inspection of the lightcurves has revealed the difference between Method I and 
Method II lightcurves: the Method II lightcurves are synchronised, while the Method 

I lightcurves are not. Lcmath processed the latter incorrectly, resulting in the artificial 
variability. 

Many PDS from combined lightcurves obtained using Method I exhibit Type A 
variability, while those obtained with Method II do not. In order to investigate the 
cause of the Type A variability in the observed PDS, we examined the Method I and 

II lightcurves of XMMU J004208. 9+4 12048 (r3-60). The Method I and Method II PDS 
are shown in the left and right panels respectively of Fig. 01 they were obtained using the 
FTOOL powspec. The Method I PDS is well modelled by a broken power law (^ 2 /dof = 
6/8). Such a PDS is characteristic of Galactic LMXBs. However, the PDS of the Method 
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FIGURE 3. Power density spectra of 0.3-10 keV lightcurves of r3-60, using Method I (left) and 
Method II (right). We present best fits to the Method I PDS of power law (PO) and broken power law 
models (BPO); a broken power law is clearly required. The PDS is hence reminiscent of Galactic LMXBs. 
However, the Method II PDS is flat, revealing that this variability is artificial. 



II lightcurve is consistent with Poisson noise; hence the variability in the Method I PDS 
is artificial. 

Following this work, we produced background-subtracted EPIC-pn lightcurves of 
various X-ray sources from the 2002 January XMM-Newton observation of the M31 
core, to see if the lightcurves from the same detector was synchronised. The Method I 
PDS again is of Type A. Hence, background subtraction also results in artefacts when 
the lightcurves are not synchronised. 

We have investigated the artefact by adding simulated lightcurves with various inten- 
sities and different delays between lightcurves [9]. We find that the artefact is determined 
by the delays, but independent of the lightcurve intensities. Hence it is vitally important 
to synchronise all XMM-Newton observations, even if the whole observation is to be 
used. 



ULTRA-LUMINOUS X-RAY SOURCES IN NGC 4559 AND 

NGC 5408 

We found published broken power law PDS in NGC 4559 ULX7 [7] and also in a ULX 
in NGC 5408 |8]. We re-analysed XMM-Newton observations of these sources, first 
following the methods described in these papers, using non- synchronised lightcurves, 
then comparing their PDS with those from synchronised lightcurves. 

Cropper et al. tm reported a broken power from an ultra-luminous X-ray source 
(ULX) in NGC 4559 with a UV/X-ray luminosity exceeding 2.1 xlO 40 erg s -1 . Fol- 
lowing their example, we extracted MOS1, MOS2 and pn lightcurves from a source 
region with 30" radius using Methods I and II. The resulting PDS are shown in the left 
and right panels of Fig. |4| respectively. The Method I PDS exhibits Type A variability, 
while the Method II PDS does not, revealing that the PDS reported for ULX 7 by 0] is 
artificial. 

Soria et al. |8] studied a series of XMM-Newton observations of a ULX in the dwarf 
galaxy NGC 5408. A Method II PDS is presented in Fig. fitting this PDS with a 
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FIGURE 4. Power density spectra of 0.2-12 keV lightcurves of the NGC 4559 ULX7 using Method 
I (left) and Method II (right). The PDS are Leahy normalised, so that Poisson noise has a power of 2. 
While a broken power law PDS is observed in the Method I PDS, there is no excess power observed in 
the Method II PDS. Hence the broken power law PDS reported for NGC 4559 ULX7 is artificial. 



Leahy Power of 2 (Poisson noise) results in a y} '/dof of 71/1 1, hence the variability is 
significant. However, there is no statistical requirement for the claimed break in the PDS 
in the ~3 ks of good data; a longer observation would be required to identify any break. 



CONCLUSIONS 

Some of the authors have previously reported variability in XMM-Newton observations 
of M31 X-ray sources that are characterised by broken power law PDS, which are sig- 
natures of disc-accreting X-ray binaries ||4j,|3j,|5|,|g]. However, these "Type A" PDS ob- 
served in XMM-Newton observations of M31 were found to be false positives, and may 
be entirely attributed to errors introduced by lcmath when combining non- synchronised 
lightcurves. It is therefore imperative to synchronise all lightcurves; this is most easily 
achieved by following Method II. 

Extragalactic analogues of the Galactic X-ray binaries should exhibit Type A and 
Type B behaviour; Type A variability should be observed in sufficiently bright X-ray 
sources. We have shown the broken power law PDS of NGC 4559 ULX7 to be artificial 
also; however, the PDS of the NGC 5408 ULX shows that timing analysis of XMM- 
Newton observations of extra-galactic X-ray sources is viable after all. Indeed, M31 
X-ray sources have exhibited pulsations iflU flHl: bursts iflHl : periodic dipping due to 



photoelectric absorption HHEj], from a precessing disc in one case ifTTIl : lastly, branch 
movement in a trimodal colour-intensity diagram reminiscent of Galactic Z-sources H 1 Oil - 
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FIGURE 5. Method II PDS for the NGC 5408 ULX, averaged over 9 intervals of 64 bins. The PDS is 
Leahy normalised, and a line of constant power equivalent to the Poisson noise was fitted; the resulting 
X 2 /dof = 71/11, showing that the source is variable. 
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